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A brief review is given of the lattice QCD calculation of the hadron spectrum. The 
status of current attempts toward inclusion of dynamical up, down and strange quarks 
is summarized focusing on our own work. Recent work on the possible existence of 
pentaquark states are assessed. We touch upon the PACS-CS Project for building our 
next machine for lattice QCD, and conclude with a near-term physics and machine 
prospects. 
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1. Introduction 

The spectrum of hadrons is a fundamental entity for lattice QCD elucidation of the 
dynamics of the strong interactions. Calculating the spectrum of standard hadrons 
such as pion, nucleon, J/ip, Y, and Z3's and B's provides a basic mean for confirma- 
tion of the theory and a benchmark for the lattice calculational methods. Finding the 
spectrum of non-standard hadrons, e.g., glueballs, hybrids and multi-quark states, 
offers chances for predictions and thereby further tests of QCD. For these reason 
the hadron spectrum has been repeatedly and routinely calculated over the years 
with increasing accuracy and sophistication. 

At present, the focus of spectrum calculations and many other subjects with 
phenomenological impact centers around dynamical simulations including the sea 
quark effects of all three light quarks, up, down and strange. There are two major 
attempts, one usingthe staggered quark actiorP and the other using the Wilson- 
clover quark actiorW In addition, a third attempt with the domain-wall quark 
action, made possible with the commissioning of QCDOC with lOTflops-class ca- 
pability, are beginnin g5l. In Sec. El we discuss results from our own attempt on the 
light hadron spectrum and quark masses, and recent preliminary results on heavy 
quark systems. 

Search for possible pentaquark states has occupied a lot of experimental and 
theoretical effort. A number of lattice QCD simulations have also been made to 
see if there is signal indicating their theoretical presence. A crucial and special 
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issue with the pentaquark concerns distinguishing bound or resonance states from 
scattering states. Recent lattice studies focused on applying theoretical criteria on 
this point, and we shall review them in Sec. [31 We conclude with a brief summary 
in Sec. E| 

2. Nf = 2 + 1 lattice QCD simulations with the Wilson-clover 
quark action 

There are two major lattice QCD research collaborations in the Tsukuba area in 
Japan, the CP-PACS Collaboration based at University of Tsukuba and the JLQCD 
Collaboration based at KEK. Since 2001 the two collaborations have jointly pur- 
sued a project to carry out simulations including dynamical up, down and strange 
quarks. The strategy adopted, and the necessary preparations carried out prior to 
actual runs, is as follows: (i) use Iwasaki RG-improvcd gluon action to span a range 
of lattice spacing toward coarse lattices and avoid the artificial phase transition 
observed for the plaquette gluon actionSI (ii) use Wilson-clover quark action with 
a fully 0(a)-improved clover coefficient calculated with the Shcrodinger functional 
methods for three dynamical flavor and (iii) apply the polynomial HMC to handle 
the strange quarlP in addition to the standard HMC for the up and down quarks 
which are treated as degenerate. 

In the table below we list the parameters of our simulations. Runs are made 
at three values of lattice spacing equally spaced in a 2 . At each lattice spacing, 
five values are taken for the degenerate up and down quark mass in the range 
fnpsl m v ~ 0.6 — 0.8, and two values for the strange quark mass at mps/mv ~ 0.7. 
Hadron masses calculated for these quark masses are fitted with a general quadratic 
polynomial of VWI quark masses m VWI — (K^ 1 — K~ 1 )/2, and are extrapolated to 
the physical point defined either by the experimental it, p and K meson masses (K- 
input) or n, p and <j> meson masses (0- input). The agreement of the lattice spacing 
determined by the two types of input as seen in Table |2] provides a simple but 
important check on the internal consistency of the Nf = 2 + 1 calculation, which is 
not realized in quenched and Nf = 2 simulations. 

2.1. Light quark sector 

In Fig. ^ we plot the continuum extrapolation of meson masses in 
quenched(triangles), Nf = 2(squares) and the present Nf = 2 + l(circles) 



Parameters of our simulations 



P 


size 


a [fm] (i^T-input) 


a [fm] (</>-input) 


trajectory 


1.83 


16 3 x 32 


0.1222(17) 


0.1233(30) 


7000 - 8600 


1.90 


20 3 x 40 


0.0993(19) 


0.0995(19) 


5000 - 9200 


2.05 


28 3 x 56 


0.0758(48) 


0.0755(48) 


3000 - 4000 
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Fig. 1. Continuum extrapolation of meson masses, compared with those for quenched and Nf = 2 
QCDI3 Note that the quenched and Nf = 2 simulations are made with the one-loop perturbatively 
0(a)-improved clover action. Thus extrapolations are made linearly in a. 



simulations^. The solid lines for the Nf = 2 + 1 data are pure quadratic fits to 
the results for the two coarse lattice spacings for which runs and measurements 
have been completed. The agreement with experiment in the continuum limit is 
encouraging, but we need further data at the finest lattice spacing for a solid con- 
clusion. 

Figure [5] shows the continuum extrapolation for the light quark masse Values 
sizably smaller than the quenched estimate as has been strongly suggested in the 
previous Nf = 2 simulations are confirmed. Again, we wait completion of the anal- 
yses at the finest lattice spacing to quote the final result for light quark masses in 
the continuum limit. 

An important issue with the analyses in the light quark sector concerns chiral 
extrapolation. Since Wilson-clover quark action involves explicit chiral symmetry 
breaking, the chiral behavior of physical quantities deviate from that in the contin- 
uum. The Wilson chiral perturbation theory® provides the procedure to work out 
the chiral behavior for finite lattice spacings for Wilson-type actions, and work is 
under way to calculate the consequence of the approach for the spectroscopic quan- 
tities including pseudo scalar meson masses, decay constants, and quark masses on 
the one hand, and to analyze data based on those results. 

2.2. Heavy quark sector 

Calculating the physical quantities of hadrons containing heavy quarks is important 
to constrain the parameters of the Standard Model and to help explore physics at 
finer scales. A serious obstacle is large systematic errors of form 0(m,Ha) where 
ma denotes the heavy quark mass, which are large at currently accessible lattice 
spacings. A standard lore for overcoming this problem is to resort to heavy quark 
effective theory approach such as NRQCD. We wish to pursue a different approach 
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Fig. 2. Continuum extrapolations of the up, down and strange quark masses obtained with the 
X-input. The data at the finest lattice is not included in the continuum extrapolations. For com- 
parison, results for quenched and Nt = 2 QCD are overlaid. 
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Fig. 3. Decay constant for (a) D a and (b) B s in q uenched C CD from relativistic heavy quark 
approaclWJ as compared to those of effective theorieJ -^ l ^ l ^ l 



in which parameters of the quark action are tuned demanding Lorentz invariance 
of on-shell quantities for arbitrary large run a to 0(aPl. 

This approach allows the 
continuum extrapolation to be taken in contrast to effective theory approaches. 

We have been testing this approach for charm and bottom systems on Nf = 2 
and quenched configurations prior to applying it to Nf = 2 + 1 QCD. In Fig . IB] 
we show results for the decay constant of D s and B s mesons in quenc hed QCEM-iJ 
and compare them with those of effective theory a 

pproacheJEGHll Mild cutoff 
dependence of the present results and consistency at finite lattice spacings indicate 
success of our approach. 



2.3. The next step 

At present the largest limitation in our data is a rather large value of the up and 
down quark mass m u d- In terms of the physical strange quark mass m p a hys , it only 
goes down to m u d ~ m p s hys /2 while experimentally m u d ks m p s hvs /25. We wish to 
reduce the value to at least m u d ~ m phvs /5 or below in order to control the chiral 
behavior. 

This requires both an enhancement of the computing power and an improvement 
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of the algorithm. The situations in both respects are very prom ising at present. For 
the latter, the domain-decomposition acceleration of HM(J^ offers a very promis- 
ing resolution. On the latter, in December 2004, Japanese Government formally 
approved our proposal to develop a massively parallel cluster at the Center for 
Computational Physics, University of Tsukuba. We list in table below the specifi- 
cation of the cluster, which is named PACS-CS (Parallel Array Computer System 
for Computational Sciences). The development is well under way, and the installa- 
tion and start of operation are scheduled in June- July of 2006. We plan to combine 
these developments to break through the current limitations in our Wilson-clover 
program. 

3. Lattice pentaquark search 

A standard lattice QCD calculation of the mass of a hadron starts with a preparation 
of the operators having the desired quantum numbers and an examination of the 
large-time behavior of the two-point function to extract the ground state signal. 
With lattice pentaquark searches, we have to distinguish a possible pentaquark 
bound state or resonance from the nucleon-kaon scattering state. Hence, we need 
to disentangle at least two states in the correlation function, and we also have to 
distinguish a bound state/resonance from scattering states. In addition, since the 
spin-parity of the pentaquark state is yet unknown, we have to explore over a large 
operator space. 

The initial lattice studies^ did not explore these points in detail. Subsequent 
calculation ^ 1 ^ 1 addressed them using two techniques. For multi-state analyses, 
a variational method using a set of operators and diagonalizing the normalized 
correlator matrix has been known for a long time. This method has been employed 
in a number of recent pentaquark searches, and has been reasonably successful. 
For distinguishing a scattering state, a basic strategy is to examine the spatial 
size dependence of the energy eigenvalues obtained in the multi-state analysis. A 
scattering state, if the interactions between the scattering hadrons are weak as in the 
case of the nucleon and kaon for the pentaquark case, would show a size dependence 



Design specification of PACS-CS 



Number of nodes 


2560 


Node configuration 


Peak performance 


14.3 Tflops 


single LV Xeon 2.8GHz 5.6Gflops 


Total memory 


5 TByte 


2 GByte memory 6.4GByte/s 


Total disk space 


0.41 PByte 


160 GBytex2(RAIDl) local disk 


Interconnect 


16 x 16 x 10 


Interconnect bandwidth 




3-dim hyper-crossbar 


250MByte/sx2/link 


OS 


Linux and SCore 


750MByte/sx2/node 


Programming 


Fortran90, C, C++, MPI 




System size 


59 racks 


Estimated power 545 kW 
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Fig. 4. Spatial size dependence of energies in the pentaquark channel in quenched qceGHI 



expected from the sum of two energies ^Jm 2 + (2tt£/L) 2 with I — 1, 2, 3, • ■ •. Another 
criteria is to look at the residue of the state in question in the two point function : 
< O(t)O(0) >^5° Zexp(-mt). For a scattering state we expect the overlap of the 
two particles decreasing as Z oc 1/V. 

In Fig. 01 we reproduce the spatial size dependence of the energies in the 1/2— 
and 1/2+ channel calculated with five operator basis in quenched QCE 18 . The 
dashed lines indicate the size dependence expected for scattering states, and the 
horizontal dotted lines show the experimental value quoted for the 6+ state. As the 
figure indicates, the authors conclude that there is no evidence for the existence of 
pentaquark states in their data. 

While conclusions are not unanimous among the studies^^l, it is our view 
that the data published to data are more consistent with the absence of pentaquark 
states at present. Caution should be stated that all data so far have been generated 
in quenched QCD, mostly at a fixed lattice spacing and at relatively large quark 
masses. Truly realistic tests with dynamical light quarks with large volume and 
small lattice spacing require further work in the future. 



4. Conclusions 

Currently lattice QCD simulations are in a state of transition. Realistic simulations 
including dynamical up, down and strange quarks are becoming routine, and chirally 
invariant quark actions are beginning to be used increasingly frequently. This trend 
will accelerate as computers with 10 Tflops-class capability, starting with QCDOC 
and followed by ApeNEXT, PACS-CS and also BlueGene/L at several institutions, 
becomes available. We hope that single hadron properties will be fully understood, 
inviting us to venture into the World beyond, in the near future. 
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